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The displacement amplitudes and the phase angles of vertically vibrating 
rigid circular plates on an elastic isotropic homogeneous half-space have 
been expressed in terms of the mass of the plate, the static spring constant 
multiplied by a frequency-dependent function, and a damping term. The 
results have been modified to apply to vibrating soil-foundation systems. 
The effects of hysteresis damping, nonlinear load-deflection characteristics 
of soils, the static prepressure, the change of soil properties with depth, and 
the difference between static and dynamic stress-strain relations of soils 
have been considered. The mathematical model has been compared with 
data on vibrations of circular foundations. The agreement between the model 
and the experimental data for cohesive soils is very good. 

I. INTRODUCTION 

Numerous attempts have been made to develop a mathematical 
model capable of representing the steady-state vibrations of a soil- 
foundation system. E. Reissner 1 solved the problem of vertical vibrations 
of a rigid circular plate on a semi-infinite elastic solid. A sign error in 
his work was discovered by O. J. Sechter 2 who presented also a corrected 
analytical solution for this case. T. Y. Sung 3 continued this work for 
different pressure distributions between the plate and the solid. G. N. 
Bycroft 4 presented approximate solutions of the steady-state vibrations 
for the degrees of freedom of a rigid circular plate on an elastic isotropic 
half -space and on an elastic stratum. Since the mathematical solutions 
become rather difficult, this approach has been used only for a strongly 
idealized soil, namely the semi-infinite elastic isotropic solid. Another 
approach to this problem is the determination of a simplified mathe- 
matical model capable of describing the vibrations of a soil-foundation 
system. O. J. Sechter 2 showed that the amplitude-frequency response 
curve of a vibrating system consisting of constant mass, viscous damp- 
ing force, and linear spring constant differs only slightly from that 
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caused by vertical vibrations of a circular plate on an elastic half-space, 
provided the mass of the plate is multiplied by a factor and the damping 
coefficient is chosen appropriately. This approach has been extended 
by M. Novak 50 who introduced a frequency independent nonlinear 
restoring force. G. Ehlers 7 attempted to simulate the soil-founda- 
tion system by assuming the foundation to be supported by a trun- 
cated pyramid of soil. This approach was used also by A. Pauw, 8 
and H. A. Balakrishna Rao and C. N. Nagaraj. 9 However, none of the 
above mentioned models was capable of representing the amplitude- 
frequency response curves of several soil-foundation systems investi- 
gated by R. K. Bernhard and J. Finelli, 10 who concluded that an equiva- 
lent system could not be analogous to a simple spring-mass system 
with viscous damping and linear elasticity. 

The purpose of this paper is to present a quasilinear system sufficiently 
adaptable to represent the amplitude-frequency response and phase 
angle-frequency curves of circular foundations on soil. 

II. THE VIBRATIONS OF CIRCULAR FOUNDATIONS ON AN ELASTIC HALF- 
SPACE* 

The vertical steady-state vibrations of circular foundations on a 
homogeneous elastic isotropic half -space have been investigated analyti- 
cally by E. Reissner, 1 0. J. Sechter, 2 T. Y. Sung, 3 and G. N. Bycroft. 4 
According to their work the displacement amplitude becomes 

„ Fo i / *' + /g2 (1) 

yo ~GFoV (l + W/i) 2 + (Mo 2 / 2 ) 2 " 

The phase angle is 

tan * --<•, + ww +/«• (2) 

where 

yo = displacement amplitude of foundation, 

Fo = force amplitude of harmonic exciting force, 

G = shear modulus of half-space, 

r = radius of circular foundation, 

b = mo/pro 3 = mass ratio, 

mo = mass of circular foundation, 

p = mass density of half-space, 

a = rotay/p/G, = frequency factor, 

* The letter symbols used in this paper are defined where they first appear and 
are arranged alphabetically, for convenience or reference, in the Appendix. 
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* = phase angle between exciting force and displacement, 
G, = dynamic shear modulus of half-space, 
w = angular frequency, 
/i and fi are functions of the dimensionless frequency factor oq , of 
Poisson's ratio v of the half -space and of the assumed pressure dis- 
tribution between the half-space and the foundation. 

T. Y. Sung 3 determined the functions f x and fz for different contact 
pressure distributions and different Poisson's ratios in terms of power 
series of the frequency factor a for a ^ 1.5. The values of these func- 
tions for a pressure distribution caused by a rigid circular plate are 
shown in Fig. 1. The functions /i and/2 as determined by G. N. Bycroft 4 
differ somewhat from those used in this paper which are based on T. Y. 
Sung's work. It should be noted that for zero frequency, /1 approaches 
a finite value and fi becomes equal to zero. 

For zero frequency, (1 ) becomes 



y . = - ^rU 



(3) 




0.6 0.8 1.0 

FREQUENCY FACTOR, a 



Fig. 1 — Functions /i and f 2 for rigid circular plates. 
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where 

y 0B = lim yo = displacement of foundation due to static load Fo , and 

6)-»0 

/io = lim /i = value of /i for zero frequency. 
<j -*o 

From Fig. 1 it follows that the value of /io can be expressed as 

-/» = if-'. (4) 

Substituting (4) in (3) results in an expression for the displacement of 
a rigid circular plate on a homogeneous elastic half -space due to a vertical 
static force. This expression is identical to that derived by Boussinesq 11 
in 1885. 

It is readily verified that (1) and (2) can be written as 

V ° V[C a a - moo; 2 ] 2 + [CSV 

tan $ = -pj 5 (6) 

C s a — Wow 2 



where 



C._-*>-J2i--fim^ (7) 



/ 10 1 — v o.-o yo 
/i/io 1 - 



/i 2 + / 2 2 4 /x 2 + / 2 

_ /2/10 1 — v fi 

/1 2 + H 4 /? + / 2 2 ' 



(8) 
(9) 



C represents the static spring constant of a rigid circular plate on an 
elastic isotropic homogeneous half-space. Its magnitude depends on 
the shear modulus and Poisson's ratio of the half-space and on the radius 
of the rigid plate. 

a and are functions of the frequency factor a and of Poisson's ratio 
v as shown in Fig. 2. They have been computed using the values of / x 
and fi shown in Fig. 1. The values of the functions a and /3 for v = and 
v = I are plotted in Fig. 2 only, because the valu es for v = \ and v = \ 
fall between those shown. Since a = row y/p/G, , the functions a and 
depend on the magnitude of the radius of the plate, r , the frequency, w, 
the dynamic shear modulus, G„, of the half-space, and the mass-density 
p, of the half -space. 
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Fig. 2 — Functions a and /3 vs frequency factor a . 



The square root of the ratio of the dynamic shear modulus, G s , to the 
mass-density, p, represents the shear wave velocity, v s , 



,V!- 



(10) 



It should be mentioned that, for many materials, the dynamic shear 
modulus, G 8 , as obtained from wave velocity measurements, is larger 
than the shear modulus, G, determined by means of static or slow dy- 
namic tests. 

For large shear wave velocities, the function a approaches one and /3 
approaches zero. For this case, (5) represents the displacement ampli- 
tude and (6 ) the phase angle of a linear vibrating system without damp- 
ing. This system consists of the mass of the plate, a weightless spring 
with a spring constant defined by (7), and a harmonic exciting force. 
It becomes apparent that the functions a and /3 modify this simple 
system to account for the additional effect of the mass-density of the 
half-space. 

The functions a and /3, as shown in Fig. 2, can be expressed with 
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sufficient accuracy as follows: 



a = 1 — eao 2 = 1 — e ( — 



@V (id 

where e and 61 are appropriately selected constants. The values of e and 
bi depend on Poisson's ratio; e varies between 0.26 and 0.30, and 
between 0.76 and 0.85. 

III. THE PARTICIPATING MASS OF THE HALF-SPACE 

H. Lorenz 12 introduced the concept of the participating soil mass in 
order to explain some results of foundation vibration tests. O. J. Sechter 2 
proved that the vertical vibration of a circular plate on an elastic half- 
space can be expressed, with sufficient accuracy, by means of a constant 
mass, a viscous damping force, and a spring constant. The mass consists 
of the mass of the plate, m , multiplied by a factor which accounts for 
the mass of the half-space vibrating in-phase with the plate. The dis- 
placement amplitudes and the phase angles of such a system can be 
written as 

Fo (13) 



yo V[c s - mrf? + IM 2 

tan<P= b "_ , (14) 

where 

m = mo + m, = constant mass, 

m, = participating mass of half-space, and 

b = damping constant. 

Equations (13) and (14) are identical to (5) and (6), provided the 
constant mass and the damping constant are selected as follows: 



m = 



— [ 1 + '=©1 

6 = &x-°C.. (16) 

v, 

This identity is readily shown by substitution of (15) and (16) in (13) 
and (14). 
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IV. DYNAMIC STIFFNESS 

Some investigators, such as W. Heukelom, 13 14 expressed the results 
of foundation vibration tests in terms of the dynamic stiffness and the 
phase angle. The dynamic stiffness is defined as the ratio of the force 
amplitude between the foundation and the subsoil to the resulting dis- 
placement amplitudes. The steady-state vibrations of a rigid circular 
plate on an elastic half-space shall now be expressed in these terms. 

The force acting at the surface of the half -space can be expressed as 

F = F - m4 (17) 

where 

p = F Q e' iwl+ * ) = force acting between plate and half-space, 

/''o = force amplitude acting between plate and half-space, 

^ = phase angle between F and the resulting displacement, 

F = Foe'' 1 "'**' = exciting force acting on plate, 

y = yoe " = displacement of plate, 

.'/ = —yawe*" = acceleration of plate, 

* = phase angle between F and the resulting displacement, 

i =V^T. 

The force amplitude acting between the plate and the half-space becomes 



F\ = VFo 2 + 2F„m i/ow 2 cos * + (m 2/o« 2 ) 2 (18) 

and for the phase angle between the force acting at the surface of the 
half-space and the displacement, the following expression is obtained: 

F Q sin <t> , in , 

tan * = -= — : j . (19) 

F cos <J> + wioJ/ow- 

The dynamic stiffness is defined as the ratio of the force amplitude 
acting between the plate and the half-space and the resulting displace- 
ment amplitude and becomes, by means of (18), 



a - * = *? J\ + 2 nfs j cos * + (ap*Y (20) 

where 

S = dynamic stiffness. 

Substituting y and <t>, expressed by means of (5) and (6), in (20), 
the dynamic stiffness becomes 

S = C, V« 2 + 2 (21) 

and by substitution of (5) and (6) in (19), the following expression is 
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obtained for the phase angle SI>: 

tan¥ = -. (22) 



a 



Substitution of the values for a and from Fig. 2 or of the approximate 
expression for a and given by (11) and (12) in (21) shows that the 
dynamic stiffness of the half-space is essentially equal to the static 
stiffness or the static spring constant, C s . The effect of a and /? on the 
dynamic stiffness is negligibly small. 

Substitution of a and /3 in (22) shows that the phase angle between 
the force acting at the surface of the half-space and the displacement 
remains always smaller than 90 degrees for values of the frequency 
factor, a , smaller or equal to 1.5. 

V. DAMPING 

Energy propagation to infinity provides the only damping of the 
vibrations of a rigid circular plate on an ideally elastic homogeneous 
half-space. The energy dissipated per cycle becomes 

D. = j>F dy (23) 

where 

D 8 = dissipated energy per cycle, 
F = Fo sin (<at + $), 
y =?/o sin tit. 

By substitution of (5) and (6) in (23) the following is obtained: 

D. = vCPyo 2 . (24) 

The energy required for the deformation of the half-space becomes 

U s = \Coty? (25) 

where 

Ut = energy required for the deformation of the half -space. 

It should be noted that the energy required for the deformation of the 
half-space is equal to the maximum energy stored in the half-space 
reduced by the inertia effects of the half-space. 

The loss coefficient of the circular plate on the ideally elastic half- 
space shall be denned as 

D. p (c?p , 

V. = ttjT = " (20) 

ZttU, a 
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where 

i) a = loss coefficient of vibrating circular plate. 

The loss coefficient of the circular plate on the elastic half-space should 
be identical to the phase angle between the force amplitude acting 
between the plate and the half-space and the resulting displacement 
amplitude shown in (22). 

It should be pointed out that both the dissipated energy and the 
energy of the half-space are proportional to the static spring constant 
ft . The loss coefficient of the vibrating system, however, is a function 
of a and /3 alone, a and (3 are defined by means of (8) and (9). Equations 
(11) and (12) are approximate expressions of a and /3. 

The medium of the half-space was considered to be ideally elastic and 
hence, the propagation of energy to infinity provided the only source of 
energy dissipation. However, engineering materials and particularly 
soils are not ideally elastic but dissipate energy due to internal friction 
if subjected to cyclic stresses. The author and R. R. Hart' 5 determined 
the specific damping capacity or the coefficient of energy absorption of 
some granular soils by means of laboratory tests. The specific damping 
capacity is defined as the amount of energy absorbed by a unit volume 
of the material per unit of energy spent for deformation per cycle. The 
specific damping capacity for these granular soils varied between 0.4 
and 0.9. D. D. Barkan 16 reported values between 0.64 and 0.79 for sand. 
For clayey soils the damping capacity varied between 0.2 and 0.6. 
Furthermore, it was established that the specific damping capacity of 
sandy and clayey soils does not depend on the rate of stress application, 
the frequency of changes in the stress, the maximum alternating stress, 
or the static prestresses. The loss coefficient is obtained by division of 
the specific damping capacity by 2ir. 

-I w 

where 

rj = loss coefficient of material 

^o = specific damping capacity of material. 

The loss coefficient, 77, is considered to be a basic property of a material. 
The difference between the loss coefficient of a material, 77, and the loss 
coefficient of a vibrating system, r) 8 , should be noted. 

G. N. Bycroft 4 investigated the effect of internal friction of the me- 
dium of the half-space on the functions /1 and j% shown in Fig. 1. /1 be- 
comes numerically smaller and / 2 greater for small values of the fre- 
quency factor, a . Changes of /i and f 2 affect the function a and shown 
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in Fig. 2 and denned by means of (8) and (9). The effect of such a change 
on the function a is rather small and a can still be expressed, with 
sufficient accuracy, by means of (11). The function j8 becomes approxi- 
mately 

= v + ha o = U + fci^. (28) 

For small values of the frequency factor, a„ , or for high shear wave veloc- 
ities and small radii of the plates, the effect of the loss coefficient of the 
medium becomes significant. 

An approximate expression of the loss coefficient of the vibrating 
circular plate is obtained by substitution of (11) and (28) in (20) 

t " + <* (-) 

V + Oitto \ v * I (on\ 



r W 



For zero frequency, the loss coefficient of the plate becomes equal to 
the loss coefficient, of the medium of the half-space. This is true only 
if the loss coefficient of the material is independent of the magnitude of 
the applied static and cyclic stresses. The author 17 measured the loss 
coefficient of circular plates on clayey silt subjected to slowly alternating 
loads. The results varied between 0.1 and 0.15. Considering the additional 
friction between the plate and the soil, the results agree well with those 
obtained by laboratory tests of sandy and clayey soils as reported by 
the author 15 and by D. D. Barkan. 16 

VI. NONLINEAR RESTORING FORCE 

The static spring constant of a rigid circular plate on an elastic iso- 
tropic homogeneous half -space has been expressed by means of (7). 
For this idealized medium, the deflection is a single-valued linear func- 
tion of the applied force. For engineering materials, and particularly for 
soils, this relationship becomes considerably more complex. Due to 
the effect of hysteresis damping, the deformation caused by the appli- 
cation of an alternating force is no longer a single-valued function. Two 
different forces cause the same deformation under the condition of 
loading and unloading. Furthermore, the average of these two forces 
may be a nonlinear function of the deformation and/or the rate of 
deformation. Other factors, such as static preloads and change of the 
properties of the half-space with depth, further complicate the force- 
deformation relation. 
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In view of the mathematical complexity of strongly nonlinear vibra- 
tion problems, it shall be assumed that the nonlinearities of the force- 
displacement function within a limited range are small. For this case 
the nonlinear force-displacement function can be expressed with suffi- 
cient accuracy by means of an equivalent linear function. 

F. = f(y) « CMv for -y Q < y ^ y<t (30) 

where 

F' = f(y) = nonlinear force-displacement function, and 
C»(yo) — amplitude dependent spring constant. 

The spring constant, C, , should be selected in such a manner that the 
stored energy at the displacement amplitude of both the nonlinear and 
the equivalent linear systems are equal. 

U = / F.dy = / f(y)dy = \C.(yn)y^ (31) 

Jo J o 

where 

U = stored energy. 

The spring constant as a function of the displacement amplitude is 
readily obtained by solving (31) with respect to C, . The static spring 
constant can frequently be expressed by the following empirical equation: 

C. = Co2/o~" (32) 

where 

Co and n are appropriately selected constants. 

It should be noted that (32) is an empirical function which has been 
used successfully by the author 17 to express the amplitude, dependent 
load-deflection characteristics of circular plates on some soils. Other 
functions may be more suitable to define these characteristics for 
different soils. 

The amplitude-frequency response curves and the phase angle fre- 
quency curves of circular plates on a quasielastic homogeneous isotropic 
half-space with hysteresis damping are obtained by substitution of 
(11), (28), and (32) in (5) and (6). For the amplitude-frequency response 
curves, biquadratic equations of the frequency as a function of the ampli- 
tudes are obtained. 

VII. THE VIBRATING SOIL-FOUNDATION SYSTEM 

The load-deflection characteristics of circular rigid plates on soil are 
generally expressed in terms of the coefficient of subgrade reaction 
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which is defined as the ratio of the average applied stress to the resulting 
displacement. The spring constant of a rigid circular plate on soil 
becomes then 

C s = irr \ (33) 

where 

k 8 = coefficient of subgrade reaction. 

By means of (7) and (33) the coefficient of subgrade reaction of cir- 
cular rigid plates on an elastic isotropic homogeneous half-space can be 
expressed as 

fc - ,, *° , ■ (34) 

7r(i — v)n 

For this case, the coefficient of subgrade reaction becomes inversely 
proportional to the radius of the plate. 

The author 17 determined the coefficient of subgrade reaction of a silty 
clay experimentally. Circular plates of different radii, r , were loaded to 
produce various average prepressures, p, between the plate and the sub- 
soil. This static prepressure was maintained until the creep settlements 
of the plates became negligibly small. An average stress, <x„ , was then 
repeatedly added to and subtracted from the prepressure, p. The result- 
ing displacements were measured. Fig. 3 shows the results schemati- 
cally. The coefficient of subgrade reaction, k„ , is calculated by dividing 
the alternating stress amplitude, a n , by the corresponding vertical dis- 
placement, y n . The results may be summarized as follows : 

(i) The coefficient of subgrade reaction decreased with an increase 

of the alternating stress amplitude, a n ■ 
(ii) An increase of the static prepressure, p, caused an increase of the 

coefficient of subgrade reaction. 
(Hi) An increase of the radius, r , caused a decrease of the coefficient 

of subgrade reaction; however, this decrease was generally smaller 

than (34 ) , valid for the elastic half -space, would indicate. 
The author and S. R. White 18 have shown that the dependence of the 
coefficient of subgrade reaction on the static prepressure can be rather 
significant for some soils. 

S. D. Wilson and E. A. Sibley 19 have shown that the shear modulus 
increases with an increase of the rate of load application. A modulus 
value determined from velocity measurements was almost an order of 
magnitude higher than that determined by means of static compression 
tests. Typical data reported by R. V. Whitman 20 show a similar phe- 
nomenon. Only for sands, the shear modulii obtained by seismic tests 
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: ALTERNATING AVERAGE 
STRESS AMPLITUDE 




VERTICAL DISPLACEMENT 

Fig. 3 — Coefficient of subgrade reaction. 



did not differ much from that determined by means of static compression 
tests. Based on these results, it appears reasonable to assume that the 
coefficient of subgrade reaction of some soils increases also with an in- 
creased frequency of the load application. 

Because the static prepressure, p, and the radius, r , of the foundation 
are important parameters affecting the magnitude of the coefficient of 
subgrade reaction, k a , the mass, m , of the foundation should be expressed 
as 



1 2. 

?Wo = - TT V 

9 



(35) 



where 

g = acceleration of gravity. 

Substituting (11), (28), (33), and (35) in (5) and (6), the following 
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expressions for the displacement amplitude, y Q , and the phase angle, $, are 
obtained: 






rrih /j/Ti - {1 + ed„ 2 } ^- « 2 1 + U + &ip«J 



(36) 



_i_ j, r ° 

tan * = - : (37) 

1 - jl + ea 2 } rP-a 



fog 



where 

ro . /fog 
Vs V v 

The resonance frequency becomes 



(38) 



fog / k>9 



. / K s g / *>y 

"" ~~ V W+7^i = . / J, , -ro'k,g-[ ■ (39) 



In order to estimate the magnitude of a , (10) and (34) are sub- 
stituted in (38) to yield 

where 
y = pg = density of soil. 

For most soils, with the exception of sand, the shear modulus obtained 
by means of slow vibration tests is smaller than that obtained by seismic 
velocity measurements and hence G/G, < 1. The specific weight of 
soils, 7, is generally smaller than the average contact pressure between the 
foundation and the subsoil, p, and therefore y/p < 1. Since the founda- 
tion is considered to be rigid, the magnitude of the radius, r , must be 
limited. An increase of the radius r would require a corresponding in- 
crease of the average contact pressure, p. Hence, it is reasonable to 
assume that d is, for many materials, a small quantity and that ea is 
then small compared to one. For this case, (39) can be simplified to 

wo « Vk.g/p. (41) 

The amplitude at resonance, y 0r , for a constant exciting force, F , fol- 
lows from (36) 
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U (42) 



yor = • r 

TTVo k» 77 + 6l — OJO 

L v* 

and for a centrifugal exciting force F = mie u 2 where m^ is the eccentric 
mass and e the eccentricity of the rotating weights: 



y° r = T Z — 1 = ~7T^ 



/r wo\ 2 ~l 



(43) 



The loss coefficient at resonance can be expressed by means of (29) and 
(43) as 






\ p, / mie g mie , 

^ = 7^ w. = — — = • (44) 

_ (r&>o\ irr pyor rmyor 

It is now hypothesized that (136) and (37) describe the amplitude- 
frequency response and the phase angle-frequency relation of a rigid 
circular foundation on soil, provided the coefficient of subgrade reac- 
tion, k s , the shear wave velocity, v, , and the loss coefficient of the soil 
y, are determined experimentally and the constants e and h\ are selected 
according to the model of the rigid circular plate on the elastic half- 
space as 0.26 and 0.7(5, respectively. For sand, with its relatively low shear 
wave velocity, the term r u)o/v a may exceed 1.5 and (36) to (39) are no 
longer applicable. 

The validity of this hypothesis can only be established by the evalua- 
tion of the data obtained by means of vibration tests. 

VIII. THE EVALUATION OP' VIBRATION TEST DATA 

Vertical vibration tests of rigid circular foundations on soil were 
conducted by H. Lorenz 1223 in Germany, D. D. Barkan 16 in the U.S.S.R., 
M. Novak 56 in Czechoslovakia, W. Heukelom 1314 in the Netherlands, 
F. J. Converse 22 of California Institute of Technology, and Z. B. Fry 21 
of the U. S. Army Engineer Waterways Experiment Station. R. K. 
Bernhard and J. Finelli 10 of Rutgers University used a rigid square 
foundation. 

A vertical sinusoidal exciting force, F, was applied to the foundation 
and the resulting displacements were measured. Generally, these dis- 
placements could be expressed with sufficient accuracy in terms of a 
sinusoidal function. 
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The weight of the foundation and its dimensions are considered to be 
known quantities. They can be expressed in terms of the average 
contact pressure between the foundation and the subsoil and the radius 
of the foundation. In case of square foundations, the radius of an equiva- 
lent circular foundation with the same contact area as the square 
foundation was used. 

The exciting force amplitude, F , is also considered to be a known 
function of the frequency. For all these reported tests, the exciting force 
has been produced by means of eccentrically rotating weights and 
therefore, can be expressed as F = m.ie u ■ 

The following experimental data can be obtained from a well con- 
ducted foundation vibration test: 

(i ) The displacement amplitude-frequency response curve, y Cw ) , and 

(it) The phase angle-frequency curve, $(w). 
However, the measurements of the phase angle-frequency curve is rela- 
tively difficult, and a number of investigators did not determine these 
data. In some cases only the resonance frequency was measured. 

Solving (5) and (6) with respect to C a a and C,/3 the following equa- 
tions are obtained: 



„ F r^oWoco . , 

C,a = — — f— + cos <S> 

Vo L Fo 



= irnfaa (45) 



CsP = -sin* = irr 2 k.p. (46) 

2/o 

k,a shall be defined as the "dynamic coefficient of subgrade reaction." 
By means of (45) and (46) it is always possible to calculate the param- 
eters C,a and C a p, provided the mass of the foundation, m , is known and 
the force amplitude-frequency relation, F (o>), the displacement ampli- 
tude-frequency relation, y (o>), and the phase angle-frequency relation, 
<t>(u>), have been determined experimentally. 

Substitution of these values for C,a and C,/3 back in (5) and (6) will 
obviously result in a perfect fit of the calculated amplitude-frequency 
response curve and the phase angle-frequency curve with the experi- 
mental data. However, according to (33), C, is a function of the coeffi- 
cient of subgrade reaction, k s . The coefficient of subgrade reaction, fc„ 
in turn is, as discussed previously, a function of the displacement am- 
plitude, y , and of the frequency, w. Due to the large amplitude changes 
occurring close to the resonance frequency, even small nonlinearities of 
the stress-displacement relation may cause a considerable change of the 
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magnitude of the coefficient of subgrade reaction. From a single ampli- 
tude-frequency response curve and the corresponding phase angle- 
frequency curve, it is impossible to determine the effect of the displace- 
ment amplitude on the parameters C,a and C$. However, if the exciting 
force amplitude is changed a number of amplitude-frequency response 
curves and phase angle-frequency curves are obtained, and it is possible 
to evaluate points of constant frequency at different displacement 
amplitudes or of constant displacement amplitudes at different fre- 
quencies. In this manner, the effect of frequency and amplitude on the 
parameters C,a and C,fi can be determined. 

The loss coefficient, tj, , of the vibrating system is defined by means of 
(26). Substitution of (45) and (46) in (26) results in 

sin <j> 

' r ' = ' (47) 



[=r+— I 



It should be noted that for phase angles $ greater than 7r/2, the denomi- 
nator in (47) becomes a difference, and small experimental errors will 
have a considerable effect on the calculated magnitude of the loss co- 
efficient, rj s . Furthermore, (47) shows again that the loss coefficient 
does not depend on the spring constant, C, . 

If the phase angle-frequency data are not available, C,a and C a /3 can 
be calculated from the amplitude-frequency response curve for specific 
frequencies only. For low frequencies, the phase angle $ can be assumed 
to be approximately equal to zero, at high frequencies equal to r, and at 
resonance equal to ir/2. Unfortunately, for a centrifugal exciting force, 
the displacement amplitudes at low frequencies are rather small and the 
data become relatively inaccurate; at high frequencies, the exciting 
force amplitude frequently exceeds the weight of the foundations and 
no data are taken. It is always possible to select frequency dependent 
functions in such a manner that an amplitude-frequency response curve 
calculated by means of (5 ) will fit the experimental data. This is possible 
for any arbitrary amplitude dependent function, C,(yo). However, there 
is no physical significance in fitting one particular amplitude-frequency 
response curve. A mathematical model capable of describing the vibra- 
tions of a vibrating soil-foundation system should not only fit the am- 
plitude-frequency response curve and the phase angle-frequency curve 
of one particular test but should also predict the effect of a change in 
radius, r , of the foundation, of a change of the static contact pressure, p, 
and of a change of the exciting force amplitude, F . 
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IX. RESONANCE DATA OF VERTICAL STEADY-STATE VIBRATIONS OF RIGID 
FOUNDATIONS 

A direct comparison of the data of vibrating rigid foundations shown 
in the literature is complicated by the different units and nomencla- 
ture used to express the results of these tests. In order to facilitate such 
a comparison, data obtained at resonance have been collected and are 
presented in Table I, using uniform units and nomenclature. 

The authors of the publications used are listed in Column 2 of the 
table. 

The type of soil is shown in Column 3. 

The radius, r , of circular foundations is shown in Column 4 and is 
expressed in feet. In case of a square foundation, such as used by R. K. 
Bernhard and J. Finelli, 10 the radius of an equivalent circular foundation 
with the same contact area as the square foundation was used. 

The average static contact pressure, p, between the foundation and the 
subsoil is shown in Column 5 and is expressed in pounds per square foot. 

Column 6 shows the angular resonance frequency, oj , expressed in 
radians per second. Resonance occurs if the phase angle between the 
exciting force and the resulting displacement is equal to tt/2 or 90 de- 
grees. For centrifugal excitation, resonance is sometimes determined by 
the tangent to the amplitude-frequency response curve through the 
origin. This is correct only for a linear damping force; however, the 
possible error is small. Even the difference of the frequency at maximum 
amplitude and the resonance frequency remains relatively small. 

The exciting force amplitude, F , is listed in Column 7 in terms of 
pounds. H. Lorenz 23 and R. K. Bernhard and J. Finelli 10 reported some 
results for which the exciting force at resonance becomes greater than 
the weight of the foundation. These data are not included in Table I. 

Column 8 shows the displacement amplitude at resonance in inches 
multiplied by 10 3 . 

The loss coefficient at resonance, a dimensionless quantity, is listed in 
Column 9 and has been calculated by means of (44). 

The compression and shear wave velocities, y c ,and, v a , are shown in 
Column 10 and are expressed in feet per second. 

Column 11 shows the coefficient of subgrade reactions, k s , which is 
defined as the ratio of the average stress caused by a static force acting 
on a rigid circular or square plate to the resulting displacement. The 
coefficient of subgrade reaction is expressed in pounds per square foot 
per inch. 

The type of data reported by the investigators are listed in Column 12. 
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Table I shows clearly the need for more complete reporting of the 
obtained test data and of associated significant soil parameters. 

Some typical values of the compression and shear wave velocities for 
different soils have been listed by D. D. Barkan 16 and are shown in 
Table II. 

X. EXPERIMENTAL VERIFICATION OF THEORY 

Equation (29) shall now be used to predict the loss coefficient, ?/„ , 
of vibrating rigid circular foundations. The results shall be compared 
with the loss coefficient, r) a , obtained by substitution of foundation 
vibration data in (47). 

Z. B. Fry 21 of the U. S. Army Waterways Experiment Station reported 
the results of vibration tests of rigid circular foundations on a silty clay. 
Fig. 4 shows a typical amplitude-frequency response curve and the 
corresponding phase angle-frequency curve. An average shear wave 
velocity v, = 475 ft /sec was determined for this soil. 

Fig. 5 shows the loss coefficient, 77,, , obtained by substitution of the 
shear wave velocity, v„ , in (29). A loss coefficient of the soil of 77 = 0.1 
was assumed, and, based on the model of the rigid circular plate on an 
elastic half-space, the constants e = 0.26 and bi = 0.76 were used. 
Furthermore, for a number of frequencies, the loss coefficient, tj 8 , was 
calculated by substitution of the data shown in Fig. 4 in (47). The 
agreement between the predicted and experimentally determined values 
of the loss coefficient is very good. The data show that the loss coefficient 
is affected to some degree by the magnitude of the exciting force am- 
plitude. 

The loss coefficient at resonance can be calculated by means of (44). 
The resonance frequencies, too, and the displacement amplitudes at 
resonance, yo r , of foundations having different radii, r , contact pres- 
sures, p, and exciting force amplitudes, F , are listed in Table I. The 
data reported by Z. B. Fry" 1 for silty clay having a shear wave velocity 
v s = 475 ft/sec are listed from numbers H-la to H-lOd. The loss co- 
efficients, t/ 8 , at resonance are shown in Fig. 6 as a function of the prod- 
uct of radius and resonance frequency, roco . The agreement between the 
loss coefficient predicted by means of (29) and these experimental 
values is excellent. 

The resonance data obtained by M. Novak are shown in Table I 
(A-la to A-lOc). The tests were conducted on a loess loam. The loss 
coefficients calculated by means of (44) are shown in Fig. 7. Unfor- 
tunately, the shear wave velocity, v s , was not reported. Therefore, a 
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Table II — Compression and Shear Wave Velocities of Some 
Soils (from D. D. Barkan) 



Soil 


Compression wave 
velocity 
[ft/sec] 


Shear wave velocity 

[ft/sec] 


Moist clay 

Loess at natural moisture 
Dense sand and gravel 
Fine-grained sand 
Medium-grained sand 
Medium-sized gravel 


4920 
2630 
1575 
985 
1800 
2460 


500 
850 
820 
360 
525 
590 



5 O.OOB 

UJ 

O 



a 0.006 




3.5 




WEIGHT OF FOUNDATION 
W = 30970LBS 

RADIUS OF FOUNDATION: 
T = 3.65FT 
SOIL: SILTY CLAY 
EXCITING FORCE 
AMPLITUDE: [LBS] 

F Q = 1.464 W 2 

F Q =1 105 0> 2 

F Q = 0.736 W 2 

F„= 0.370 W 2 
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Fig. 4 — Typical amplitude-frequency response curve and phase angle-fre- 
quency curve of rigid circular foundation on silty clay. 
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Is 



W = 30,970 LBS 

T = 3.65 FT 
• F = m,e W 2 = 1.464 CU 2 (LBS) 

n F =m,e w 2 = i.i05 (u 2 (lbs) 
o F = m,e w 2 = o.736 <D 2 (lbs) 

A F = m,e 0> 2 = 0.370 W 2 (LBS) 



300 400 

r (i> (FT. RAD/SEC) 



Fig. 5 — Loss coefficient of a vibrating circular foundation on silty clay. 
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Fig. 6 — Loss coefficients at resonance of circular foundations on silty clay. 
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number of different shear wave velocities were substituted in (29) and 
the resulting curves are also shown in Fig. 7. Table II shows a value of 
the shear wave velocity of loess v, = 850 ft/sec. The agreement of the 
predicted values of the loss coefficient and the experimentally deter- 
mined values at resonance are again very good. 

Z. B. Fry 21 reported also the results of vibration tests on sand. The 
resonance data, the radii, and the contact pressures are listed in Table 
I (I-la to I-9d). A typical displacement amplitude-frequency response 
curve and the corresponding phase angle-frequency curve is shown in 
Fig. 8. The loss coefficients, rj a , were calculated by means of (47) and 
are shown in Fig. 9. No loss coefficient was predicted because the ex- 
perimentally determined shear wave velocity decreased with an increase 
of the frequency. Fig. 10 shows the loss coefficients at resonance of 
foundations having different radii and contact pressures. The loss co- 
efficient at resonance appears to be approximately constant. Fig. 11 
shows the loss coefficients calculated by means of (47) for two founda- 
tions having different radii. The figure shows that for this sand, the loss 
coefficient is independent of the radius of the foundation. In order to 
obtain the constant loss coefficient at resonance shown in Fig. 10, the 
resonance frequency must be constant too. The resonance frequencies 
listed in Table I for this case are indeed almost constant. 




300 400 500 

r io (ft. rad/sec) 



Fig. 7 — Loss coefficients at resonance of circular foundations on loess loam. 
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RADIUS OF FOUNDATION: 
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amplitude: [lbs] 
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Fig. 8 — Typical amplitude-frequency response curve and phase angle-fre- 
quency curve of rigid circular foundation on sand. 
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Fig. 9 — Loss coefficient of a vibrating circular foundation on sand. 
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Fig. 10 — Loss coefficients at resonance of circular foundations on sand. 
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Fig. 11 — Loss coefficients of two different foundations on sand. 



It may be concluded that the proposed mathematical model is capable 
of predicting the loss-coefficient, j] a , of a vibrating circular rigid founda- 
tion on a cohesive soil. For sand, however, this model is not applicable, 
and a different mechanism must be assumed in order to explain the 
independence of the loss coefficient from the magnitude of the radius of 
the foundation. The assumption of a dynamic arching effect causing the 
vibration of a large sand mass of a definite radius may provide an in- 
tuitive explanation. 

The resonance frequency of the proposed mathematical model can be 
determined by means of (39). The resonance frequency depends on the 
coefficient of subgrade reaction, k s , the static contact pressure, p, the 
radius of the foundation, r , and the shear wave velocity, v a . Equation 
(39) can frequently be replaced by (41), which neglects the inertia 
effects of the soil. It has been shown that these effects are relatively 
small for cohesive soils and foundations of small radii. The coefficient 
of subgrade reaction at the resonance frequency is generally slightly 
higher than the statically determined value. Hence, neglect of the effect 



A VIBRATINC SOIL-FOUNDATION SYSTEM 



209 



of the rate of loading in the numerator of (39) is often compensated by 
neglecting the inertia effects in the denominator. 

The assumption, that the resonance frequency can be calculated by 
means of (41), is supported by the data obtained by D. D. Balkan 18 
shown in Table I (B-la to F-4). The calculated values of the resonance 
frequency, based on static load-deflection tests, agree well with the ex- 
perimentally determined values. The error in percent is shown in Fig. 12 
as a function of the radius of the foundation. The error does not exceed 
20 percent. In general, the predicted values of the resonance frequency 
are higher than the experimentally determined ones. This was to be ex- 
pected because the inertia effects of the soil are neglected in (41 ) . Pos- 
sible experimental errors may have caused the four predicted resonance 
frequencies shown in Fig. 12 which were lower than the experimentally 
determined values. 

Unfortunately, to the best of the author's knowledge, only D. D. 
Barkan 1 determined the static coefficient of subgrade reaction in con- 
nection with vibration tests of circular foundations. However, the dy- 
namic coefficient of subgrade reaction can be calculated by means of 
(45) provided the exciting force, displacement amplitude, and the phase 
angle are known experimentally determined functions of the frequency. 
The data shown in Fig. 4 were substituted in (45) and the dynamic 
coefficient of subgrade reaction was calculated for three different dis- 
placement amplitudes. The results are shown in Fig. 13. There is an 
expected scattering of the data. However, it may be concluded that the 
dynamic coefficient of subgrade reaction, k a a, decreases with an increase 
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Fig. 12 — Error between calculated and measured resonance frequencies. 
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Fig. 13 — Dynamic coefficients of subgrade reaction as a function of the angu- 
lar frequency for different displacement amplitudes on silty clay. 

of the frequency, a, at a constant displacement amplitude, y . Changes 
of the displacement amplitudes result, however, in significant changes 
of the coefficient of subgrade reaction. 

Fig. 14 demonstrates the dependence of the dynamic coefficient of 
subgrade reaction on the displacement amplitude. The dynamic co- 
efficient of subgrade reaction decreases with an increase of the displace- 
ment amplitude. For a constant frequency, the coefficient of subgrade 
reaction as a function of the displacement amplitude can be expressed 
with sufficient accuracy by the following empirical equation: 



/Co 



2/o 



(48) 
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where 

n = 0.26 = experimental constant, 

h, = coefficient of subgrade reaction at displacement amplitude 

yo*, and 
tjo* = displacement amplitude used as reference amplitude. 

The maximum displacement amplitude at resonance, y 0r = 0.0182 in., is 
shown in Fig. 4 or in Table I, No. 5d and shall be used as reference 
amplitude (ijo* = yor). Resonance occurred at to * = 95.4 Rad. sec -1 . 
Solving (39) with respect to the coefficient of subgrade reaction, one 
obtains 



ICa — 



gv 



(49) 
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placement amplitudes for different frequencies on silty clay. 
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average contact pressure, p*, and the radius, r *, of the reference founda- 
tion are known quantities. 

By means of (11), (48), and (49) the dynamic coefficient of subgrade 
reaction k a a can be written as follows 



k„a = 



p*Wo* 2 



[■-•( 



1 - e 



/ VoW Y"] \ y° 



(50) 



Equation (50) is plotted in Fig. 13 for four values of the displacement 
amplitude, y , and in Fig. 14 for two values of the angular frequency, w. 
A comparison of the data obtained by means of (45) with the plot of 
(50) shows that the results are compatible; however, the scattering of 
the experimental data prevents a direct verification of (50). 

So far, the radius of the foundation, r , and the subsoil were assumed to 
remain constant. Now, however, the effect of a change of the radius 
and the contact pressure on the natural frequency of rigid circular 
foundations shall be considered. From (39) follows 
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(51) 



The asterisk indicates known quantities used for a reference founda- 
tion. 

The coefficient of subgrade reaction is, as mentioned before, a function 
of the displacement amplitude, the radius of the foundation, the contact 
pressure, and the rate of loading. The author" proposed the following 
expression for the coefficient of subgrade reaction subjected to fluctuating 
loads 






-(0"[<- "(ST +'Wi 



(52) 



where n, r, c*, and p are experimentally determined constants. If these 
constants are assumed to be zero, then (52) is proportional to (34) which 
is valid for the case of the elastic half-space, n indicates the amplitude 
dependent nonlinearity of the coefficient of subgrade reaction. An in- 
crease of the stiffness of the soil with depth is considered by the term 
(r */»'o) u ~ ,) f and the effect of the contact pressure is indicated by (p/p*) p . 
It shall be assumed that the coefficients of subgrade reaction for flue- 
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tuating loads were determined at a frequency in the range of the reso- 
nance frequencies of the vibrating foundations. Hence, no provision is 
made in (52) for the effect of the rate of loading. Substituting (52) in 
(51) one obtains 
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(53) 



Fig. 15 shows the angular resonance frequencies of the rigid circular 
foundations reported by Z. B. Fry 21 for silty clay as a function of the 
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Fig. 15 — Resonance frequencies of rigid circular foundations on silty clay. 
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displacement amplitude at resonance. These data are listed in Table I, 
No. H-la to H-lOd. Z. B. Fry 21 reports furthermore, that the shear 
modulus determined by seismic methods does not change significantly 
with depth below the surface for this silty clay. Hence, it may be as- 
sumed that the constant, r, in (52) and (53) becomes approximately 
zero. Furthermore, it can be assumed that the effect of the static contact 
pressure can be neglected, since the overburden does not increase the 
shear modulus significantly. Hence, c* is assumed to be small compared 
to unity. With these assumptions, (53) simplifies to 
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95.4 Rad/sec = angular resonance frequency of reference 
foundation for maximum exciting force listed in Table I, No. 
H-5d, 

0.0182 in. = amplitude at resonance of reference foundation, 
3.65 ft = radius of reference foundation, 
738 lb/ft 2 = average static contact pressure of reference 
foundation, 
0.26, 

386 in./sec 2 = acceleration of gravity, 
475 ft/sec = shear wave velocity, 

0.26 = experimental constant determined from Fig. 14, and 
coefficient of subgradc reaction for reference foundation calcu- 
lated by (49). 



For the contact pressure p = p* — 738 lbs/ft 2 and the radius r = 
r Q * = 3.65 ft, the angular resonance frequency of the reference foundation 
is obtained and shown in Fig. 15 as a function of the displacement 
amplitude at resonance, ?/or . This curve fits the corresponding experi- 
mental data very well, which was to be expected from the results shown 
in Fig. 14. For a contact pressure p = 612 lbs/ft 2 and a radius r = 3.65 
ft the maximum values of the angular resonance frequencies of the 
here considered foundations are obtained. For p = 612 lbs/ft 2 and r = 
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5.17 ft follow the minimum values. They are also shown in Fig. 15 as 
functions of the displacement amplitude at resonance. The experimental 
data and the curves computed by means of (54) show good agreement. 
It should be realized that the changes of the radii and the contact pres- 
sures were relatively small and their effects on the resonance frequency 
appear to be in the same order of magnitude as the possible experi- 
mental errors of measuring the resonance frequency. However, the 
dependence of the resonance frequency on the displacement amplitude 
is clearly demonstrated. 

The resonance frequencies obtained by M. Novak on a loess loam are 
listed in Table I, No. A-la to A-lOc. They are shown in Fig. 16 as a 
function of the displacement amplitude at resonance. Again it is appar- 
ent, that the resonance frequency decreases with an increase of the dis- 
placement amplitude; however, not as rapidly as in Fig. 15. Unfor- 
tunately, no detailed description of the soil properties is available; 
however, it has been shown in Fig. 7 that an assumed shear wave veloc- 
ity v a = 800 ft/sec allowed the prediction of the loss coefficient with 
good accuracy. Hence, the term v(r /v e ) 2 becomes small compared to 
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one and can be neglected in (53). Equation (53) can then be written 

3 -ter ©'[<■ -«&r+'m<» 

For each set of data shown in Fig. 10 follows n = 0.18. Assuming c* « 1 
and p = 0.38, (55) becomes 

'-*\]0— Pl / *\0.09 /-*\0.31 

(56) 



a .= (yj!-Y 2 (£T* = (v«-X M (PIT 1 

0* \y0r) \pj \y0rj \P/ 



Equation (56) shows the effect of a change of the displacement ampli- 
tude and of the contact pressure on the resonance frequency for the 
data shown in Fig. 15. The agreement between the resonance frequencies 
calculated by use of (56 ) and the experimental data is excellent. 

The assumption of c* « 1 in (56) results in a resonance frequency 
which is not affected by a change of the radius of the foundation. Since 
this is possible only if the overburden increases the stiffness of the soil, 
the stiffness of the soil must also be affected by the magnitude of the 
contact pressure p. Hence, an experimental constant p = 0.38, which 
indicates the effect of the contact pressure on the stiffness of the soil, 
is reasonable. 

The effect of a change of the radius of the foundation on the resonance 
frequency shall now be discussed. The resonance frequencies of the 
foundations described by (54) decrease with an increase of the radius. 
The decrease of the resonance frequency is approximately proportional 
to the square root of the radius. A similar relationship between the 
resonance frequency and the radius of the foundation was established 
by G. P. Tschebotarioff 24 who evaluated the performance records of a 
limited number of existing engine foundations. However, foundations, 
which are described by (56), are independent of the radius. H. Lorenz 23 
reported an increase of the resonance frequency due to an increase of 
the radius. This apparent contradiction is readily explained by the 
proven nonlinearity of the coefficient of subgrade reaction. The exciting 
force in all test results here reported was produced by eccentrically 
rotating weights. Therefore, the ratio of the displacement amplitudes at 
resonance becomes, by means of (43), 
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Substituting (57) in (52) and solving for the ratio of the displacement 
amplitudes at resonance and substituting this ratio in (53) the follow- 
ing expression for the resonance frequency is obtained: 
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Since n « 1, (ijor*/uor)", and Ci" /a are approximately equal to one. For 
c* = 0, (58) becomes 
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The resonance frequencies of vibrating foundations which are expressed 
in terms of (59) are affected by the nonlinearity of the coefficient of 
subgrade reaction. For relatively small foundation radii and relatively 
large nonlinearities expressed by the experimental constant, n, an 
increase of the radius may result in an increase of the resonance fre- 
quency. For large foundation radii and small nonlinearities, however, 
the frequency should decrease for an increase of the radius. Resonance 
frequencies expressed in terms of (60) should show an increase of the 
resonance frequency with an increase of the radius. 

Equations (59) and (60) show also the effect of the contact pressure, 
p, on the resonance frequency. Generally, the resonance frequency de- 
creases for an increase of the contact pressure, however, at a rate less 
than expected from purely elastic considerations. The nonlinearities 
expressed by the experimental constant, n, in (59) and (60) reduces the 
effect of the contact pressure which is proportional to the mass of the 
foundation. The stiffening effect of the contact pressure is expressed 
by the experimental constant, p. The value of p = 0.38 was used to 
express the data shown in Fig. 16. In the plate bearing tests conducted 
by the author and S. R. White 18 , a value of p = 1.859 was obtained. For 
such a soil, an increase of the resonance frequency with an increase of 
the static contact pressure should be expected. In other words, the 
increase of the mass of the vibrating system is more than compensated 
by the resulting increase of the stiffness. 

It may be stated, that the proposed mathematical model is not only 
not contradicted by the available experimental data obtained on silty 
clay and loess loam, but that there is a satisfactory amount of evidence 
supporting this theory. 

Equation (45) was also used to calculate the dynamic coefficient 
of subgrade reaction for sand. The amplitude-frequency response curve 
and the phase angle-frequency curve shown in Fig. 8 were used for the 
calculation. The dynamic coefficient of subgrade reaction, k 8 a, divided 
by the contact pressure, p, and multiplied by the acceleration of gravity 
is shown in Fig. 17 as a function of the frequency. The data show a strong 
dependence on the frequency. The dynamic coefficient of subgrade reac- 
tion decreases rapidly if the frequency is increased. This phenomenon 
was to be expected, because the loss coefficient increased rapidly with an 
increase of the frequency as shown in Figs. 9 and 11. However, the shape 
of the dynamic coefficient of subgrade reaction-frequency curve does 
not agree with the shape of the function a shown in Fig. 2. Hence, the 
proposed mathematical model of the vibrating rigid circular foundation 
is not applicable to foundation-sand systems. The dynamic coefficient 
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of subgrade reaction was also determined for two additional foundations. 
The amplitude-frequency response curves and the phase angle-frequency 
curves of these foundations were also provided by Z. B. Fry 21 of the 
U.S. Waterways Experiment Station. The dynamic coefficients of 
subgrade reaction divided by the contact pressure and multiplied by the 
acceleration of gravity for these additional foundations are also shown in 
Fig. 17. There appears to be a difference of the dynamic coefficient of 
subgrade reaction divided by the contact pressure for low frequencies. 
For higher frequencies, however, the data are practically identical. This 
implies a coefficient of subgrade reaction proportional to the contact 
pressure and independent of the radius of the foundation. The inter- 
section of the o) 2 -curve with the data indicates the resonance frequency. 
For the three sets of data shown in Fig. 17, the resonance frequency 
appears to be constant. In order to verify this assumption, all the reso- 
nance frequencies reported by Z. B. Fry 21 are shown in Fig. 18 as a func- 
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Fig. 17 — Dynamic coefficients of subgrade reaction as a function of the fre- 
quency on sand for three different foundations. 
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tion of the displacement amplitudes. For all practical purposes, the 
resonance frequencies for these circular foundations on sand are constant, 
independent of displacement amplitude, radius, and contact pressure. 
The data by F. J. Converse 22 listed in Table I, No. J-la to J-7s, show 
a dependence on the applied force amplitude but appear to be relatively 
independent of the contact pressure and the radius. A mechanism 
describing the vibrations of circular foundations on sand must be differ- 
ent from the here proposed mathematical model. Arching of sand has 
probably a significant effect on the vibrations of foundation-sand sys- 
tems. For engineering purposes, it appears to be sufficient to determine 
the resonance frequency and the loss coefficient with a small vibrator 
and to assume an identical resonance frequency and loss coefficient of 
the planned foundation. The suggestion by G. P. Tschebotarioff 25 
of providing cavities near the four corners of machine foundation blocks 
which could be filled, or emptied, to change the weight of the block by 
15 to 20 percent would have no or only an insignificant effect on the 
resonance frequency and the loss coefficient of foundations on sand. 

The effective mass of a vibrating rigid circular plate on an elastic 
isotropic homogeneous half-space is determined by (15). The mass of 
the half-space vibrating in phase with the plate becomes 



_ rr [ T J> 



(61) 
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where 

m 8 = mass of half-space vibrating in phase with plate. 

The ratio of the soil mass to the mass of the plate becomes, by means 
of (61), (33), and (35), 

e**fcY. (62) 

m p \Vs/ 



m. _ k s g /r 
m 

Substitution of (49) in (62) gives 



g=- y ';\, - (63) 

m o . /rotgo\ 

8 \ v, ) 

The validity of (63) is limited to values of r&)o/v, ^ 1.5. For this value, 
the ratio of the participating soil mass to the mass of the vibrating plate 
becomes 1.41. For most vibrating foundations listed in Table I, this 
ratio becomes considerably smaller and the participating soil mass, m s , 
becomes insignificantly small for small foundations. This conclusion is 
apparently contradicted by the theoretical and experimental investiga- 
tions of W. Heukelom. 1314 However, it is readily shown, that the 
amplitude dependent nonlinearity of the soil, which he recognized but 
neglected, can account for this apparent contradiction. The dynamic 
stiffness can be calculated by means of (20) from experimentally de- 
termined displacement amplitudes and phase angles. The vibration 
data obtained by Z. B. Fry 21 shall be used again to calculate the dynamic 
stiffness. The displacement amplitude-frequency curve as well as some 
values of the dynamic stiffness, calculated by means of (20), are shown 
in Fig. 19. No provision was made to account for the nonlinearity of the 
load -deflect ion characteristic of the soil. W. Heukelom 1314 fitted to such 
experimental data the following function: 

S — 8 | u =o — MiO) for & ^ coo 

= S | u= o — M a [ — ) w" for d)o ^ a). 

The constants *S | u= o and m, , representing the static stiffness and the 
participating soil mass, can be selected to fit the experimental data very 
well indeed. This equation would indicate a constant participating soil 
mass, m, , for frequencies below the resonance frequency and a partici- 
pating soil mass, m, (wo/w) 8 , decreasing with an increase of the frequency 
for frequencies greater than the resonance frequency. The magnitude of 
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table i. no: h-3c 
soil: silty clay 
radius: t = 2.58 ft 
contact pressure: p=i225lb/f 
resonance frequency: 

0) = 78.5 RAD/SEC 
AMPLITUDE AT RESONANCE'. 
Y or = 49.5 XIO"" 3 IN 
EXCITING FORCE AMPLITUDE AT 
RESONANCE : F op = 6820 LBS 


T2 
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Fig. 19 — Effect of the displacement amplitude on the dynamic stiffness. 
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the participating soil mass becomes considerably larger than that ob- 
tained by (63). 

The dynamic stiffness of the proposed mathematical model is obtained 
by substitution of (33), (28), and (11) in (21) 



'■^n'-feM'^f?)! (64) 

and by substitution of (48) as 






I ;(^" 2 



+ 



t + »i(»)j. «*) 



where k s * is determined from (49). For the resonance frequency, irro z k s 
becomes 420,000 lbs/in. The displacement amplitude at resonance 
ijor = 2/o* = 0.0495 in. and the experimental constant n = 0.26 arc used 
for the calculations. The solid line in Fig. 19 indicates the dynamic stiff- 
ness for a constant displacement amplitude y = yo* = ijor ■ As is to be 
expected, this curve fits the data only at the resonance frequency. The 
dotted line, however, shows the corrections made to account for the 
amplitude dependence of the dynamic stiffness expressed by (z/o*/z/o)"- 
The displacement amplitude, yo , at each frequency is taken from the 
displacement amplitude-frequency curve also shown in Fig. 19. The 
dotted curve fits the experimental data very well. It should be noted 
here that this curve was rationally developed and not obtained by direct 
fitting of data. It may be concluded that the reduction of the dynamic 
stiffness close to resonance is caused by the amplitude dependent non- 
linearity of the load deflection characteristic of the soil-foundation sys- 
tem and not by a frequency dependent mass or a frequency dependent 
coefficient of subgradc reaction. Even a small nonlincarity will cause a 
significant change of the dynamic stiffness close to resonance due to the 
large amplitude changes occurring in this frequency range. Again, the 
proposed mathematical model is applicable only for cohesive soils and 
no attempt was made to calculate the dynamic stiffness of sand-founda- 
tion systems. 

XI. SUMMARY AND CONCLUSIONS 

The displacement amplitudes and the phase angles of vertically 
vibrating rigid circular foundations on soil have been expressed in 
terms of the average static contact pressure between the foundation and 
the subsoil, the coefficient of subgrade reaction for fluctuating loads, the 
magnitude and frequency of the exciting force, the radius of the founda- 
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tion, the loss coefficient of the soil, the shear wave velocity, and two 
theoretically derived constants accounting for the inertia effects of the 
soil. The nonlinear stress-deflection characteristic of the coefficient of 
subgrade reaction as well as the dependence of the coefficient of sub- 
grade reaction on the radius of the foundation and the average contact 
pressure between the foundation and the subsoil have been considered. 
The derived mathematical model describes, as a limiting idealized case, 
the vibrations of rigid circular plates on an elastic isotropic homogeneous 
solid. Furthermore, it is capable of expressing the static load -deflection 
characteristics of rigid circular foundation on soil and, as a limiting 
case, on an elastic half-space. 

Based on this mathematical model, the following behavior of a vibrat- 
ing soil-foundation system should be expected: 

(i) An increase of the exciting force amplitude results in a decrease 
of the resonance frequency of the vibrating foundation. This 
effect is caused by the nonlinearity of the coefficient of subgrade 
reaction. 
(ii) An increase of the radius of the foundation when maintaining a 
constant contact pressure and exciting force can result in either 
an increase or a decrease of the resonance frequency of the 
vibrating soil-foundation system. For large foundations an 
increase of the radius will generally result in a decrease of the 
resonance frequency. For smaller foundations the nonlinearity 
of the coefficient of subgrade reaction and the increase of the 
stiffness of the soil with depth may cause an increase of the 
resonance frequency for an increase of the radius of the founda- 
tion. For small nonlinearities of the coefficient of subgrade 
reaction and small increases of the stiffness of the soil with 
depth., an increase of the radius of the foundation will cause a 
decrease of the resonance frequency, 
(m) An increase of the contact pressure will generally lower the 
resonance frequency, however, less than would be expected from 
purely elastic considerations. The nonlinearity of the coefficient 
of subgrade reaction and the increase of the coefficient of sub- 
grade reaction with an increase of the applied static contact 
pressure explain this phenomenon. For some soils an increase 
of the contact pressure could possibly result in an increase of 
the resonance frequency if the coefficient of subgrade reaction 
increases very rapidly with an increase of the static contact 
pressure. 
(iv) The loss coefficient, expressing the damping of the vibrating 



A VIBRATING SOIL-FOUNDATION SYSTEM 225 

system, increases with an increase of the radius of the foundation. 

The inertia effects of the soil or of an elastic half -space account 

for this behavior. 
(v) The displacement amplitudes increase less than proportional 

to the amplitudes of the exciting force. Again, the nonlinearity 

of the coefficient of subgrade reaction explains this phenomenon. 
The predicted behavior of a vibrating soil-foundation system is well 
supported by the available experimental evidence. An evaluation of 
experimentally determined amplitude-frequency response curves and 
phase angle-frequency curves support, furthermore, the validity of the 
proposed mathematical model for cohesive soils. It is shown, that the 
soil mass vibrating in-phase with the foundation is considerably smaller 
than was previously expected. It should be noted, however, that this 
mathematical model is not applicable to foundations on sand. For this 
case, the resonance frequency and the loss coefficient, remain practically 
constant, independent of changes of the radius of the foundation or the 
contact pressure. 

APPENDIX 

Notation 

The following letter symbols have been adopted for use in this paper: 
(', = static spring constant of rigid circular plate on elastic half-space 

or load-displacement relation of circular foundation on soil. 
Co = experimentally determined constant used to express the non- 
linearity of load-displacement relation of circular foundations on 
soil. 
D„ = dissipated energy per cycle. 

F = Fffi %{f * ' = harmonic exciting force acting on plate or foundation. 
F = amplitude of exciting force acting on plate or foundation. 
F s = f(y) = nonlinear force-displacement function. 
F = F a e ,(u " +iy) = force acting between plate or foundation and half- 
space or subsoil. 
F = force amplitude acting between plate or foundation and half- 
space or subsoil. 
G = shear modulus of elastic half-space. 
G, = dynamic shear modulus of half-space. 
S = Fo/yo = dynamic stiffness. 
U = stored energy. 
U, = energy required for the deformation of the half -space. 
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a = r w -y/p/Gs = frequency factor. 
a 0r = r Q u -y/p/Gs = frequency factor at resonance. 

do = (r /v t ) -s/Cs/mo = (r /v 8 ) -\/k a g/p. 

b = damping constant of linear vibrating system. 

bi = 0.76 = constant used to approximate function /3. 

6o = mo/ pro = mass ratio. 

C* = experimental constant used to express the coefficient of subgrade 
reaction. 

e = base of natural logarithms. 

e = eccentricity of rotating masses of vibrator. 

/i = function of the frequency factor and Poisson's ratio. 
f ]0 = lim/i 

ao-0 

/ 2 = function of the frequency factor and Poisson's ratio. 
g = acceleration of gravity. 

i = V^l- 

k a = coefficient of subgrade reaction. 
Wi = eccentric mass of vibrator. 
m = mass of circular plate or foundation. 
m = mo + m a = effective vibrating mass. 
m, = participating soil mass. 
n = experimental constant used to express the nonlinear load-dis- 
placement relation of circular foundations on soil. 
p = average static contact pressure between foundation and subsoil. 
p = experimental constant used to express the increase of the coeffi- 
cient of subgrade reaction with an increase of the static contact 
pressure. 
ro = radius of rigid circular plate or foundation. 
r == experimental constant used to express the increase of the stiffness 

of soil with depth below ground level. 
t = time. 

v s = shear wave velocity of soil. 
v c = compression wave velocity of soil. 
y = displacement of plate or foundation. 
y = displacement amplitude of plate or foundation. 
y = acceleration of plate or foundation. 
ijos = displacement of foundation due to static load. 
y 0r = displacement amplitude of plate or foundation at resonance. 
$ = phase angle between exciting force and displacement of plate or 

foundation. 
S^ = phase angle between force acting between the plate and the half- 
space and the resulting displacement. 
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^o = specific damping capacity of soil, 
a = /i/io/(/i 2 -\-fi) = function expressing the inertia effect of the 

half-space on the restoring force. 
ao = a(aor). 

6 = /2/10/ (fi + /s ) = function expressing energy losses due to wave 

propagation. 

7 = density of half-space or soil. 

p = constant used to approximate function a. 
7j = loss coefficient of soil. 

y* = D„/2irU s = loss coefficient of vibrating plate or foundation. 
v = Poisson's ratio of half-space or soil, 
p = mass density of half-space or soil, 
oj = angular or circular frequency of vibrating system, 
too = angular resonance frequency. 

* The asterisk is used as a superscript to designate fixed known or 
experimentally determined quantities. 
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